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Abstract 

An experimental study of the applicability of mechanics equations 
to describing the process of equilibrium establishing in an isolated spin 
system was performed. The time-reversion effects were used at the 
experiments. It was demonstrated, that the equations of mechanics 
do not describe the experimental results. The demonstration of the 
incomleteness of quantum mechanics description of the macrosystems 
annuls the contradiction which lays in the basis of the irreversibility 
problem. 

1 Introduction 

The problem of correlation between reversibility and irreversibility in the 
evolution of physical systems is one of the fundamental problems of contem- 
porary science. It is basically the problem of the unsolved contradiction 
between the determined world that we get if we absolutize the mechanics 
equations , and the real world as we know it. At the level of the law of 
physics, this contradiction appears in the incompatibility of the 2nd law of 
thermodynamics and the reversibility of physical systems evolution followed 
from classical and quantum mechanics equations. 
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In a numerous attempts to arrive to the irreversibility from the reversible 
equations there were an avowed or hidden assumptions that would allow the 
authors to arrive to a desirable result but would not have grounds within the 
mechanical theory. 

The reversibility of the evolution of those systems that are described by 
the mechanics equations has given birth to the idea that the source of the 
irreversibility is the interaction of the system with the surrounding environ- 
ment where there are statistical laws in operation. The first thing that comes 
to mind when one reads the works proclaiming this point of view is the ques- 
tion of where those statistical laws in the physical objects surrounding the 
analyzed system come from. 

The irreversibility problem can neither be solved within the dynamical 
chaos theory, since the chaoc in question is a determined and reversible one. 

The classical mechanics was so spectacularly successful and its predictions 
were so exact, that it didn't leave great Boltzmann the only and the decisive 
argument in the discussion with his opponents - he couldn't have made a 
statement about the classical mechanics theory not being complete. 

Until early 20th centure there were no experimental data demonstrating 
the incompleteness of the classical mechanics. The study of black body ra- 
diation, atom spectra and photoeffect has lead to the creation of quantum 
theory. But with all the empressive achievements of quantum mechanics, 
there is our common sense and the fact, that the heat always goes from 
a hot object to a cold one, which prevent us from deeming the quantum 
mechanics absolutely exact. This makes it crusual to obtain experimental 
data which would not be describable on the basis of quantum mechanics and 
would demonstrate the manifestation of the 2nd law of thermodynamics at 
the evolution of macrosystems. 

It is practically impossible to solve the mechanics equations for the macrosys- 
tems in which the thermodynamics laws apply. But if the mechanical state 
of the system cannot be described, how can we compare the predictions of 
mechanics with those of thermodynamics? Nonetheless, such comparison 
is possible, and it turns out that it does not require solving the mechanics 
equations. This comparison can be made on the basis of real experiments 
in changing the time sign in the equations which describe the macrosystem 
behaviour. 
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1.1 Possibility of time reversion experiments 

It is possible to speak about a time-reversing experiment in principle based on 
the fact that changing the system Hamiltonian sign is equivalent to changing 
the time sign. This statement becomes obvious if we present the system wave 
function as 

tf(t) = e _iW **(0) (1) 
or write down the Liouville equation solution for the system density matrix: 

a(t) = exp(-iHt)a(0) exp(iHt) (2) 

We should bear in mind, though, that we are interested in the internal 
evolution of the system, which is determind by the interactions between all 
its particles. It might seem impossible to change the sign of interactions 
in macrosystems where the thermodynamics laws apply. Nevertheless, at 
least for spin systems, it turned out to be possible not only to change the 
Hamiltonian sign, but also to compare the forecasts of mechanics and the 
effects of the 2nd law of thermodynamics under those conditions. 

The Hamiltonian of a spin system placed in a strong external constant 
magnetic field H = in a frame rotating with a frequency c^o around the 

axis Z reads ^: 

W = H' A = ^A^Jzj - + I-il+j)] ( 3 ) 

i<j 

If an alternating field of a resonant frequence and an amplitude 0^1/7 is 
applied to the system, then in the rotating frame we get: 

H R = u> 1 I x + H' d (4) 

Going into the tilted rotating frame by a unitary transformation deter- 
mined by the operator exp(— i\I y ) , we obtain: 

Htr = u x I z - \H-'d + f -P> 

p = nf + n d ' 2 \ nf = n d - 2> = e a tJ i +l i +r (5) 

Kj 
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According to (2), we write down for the density matrix of the system with 
the Hamiltonian (5): 

a{t) = expi-ifal, - \n' d + \p)}a{Q) exptfM, - hi' d + jjp)} (6) 

In a strong alternating field Hi, when uj\ ^> ool (wl = where 
the local field H L = {[Tr(H' d ) 2 ]/[Tr(M^)]}^ 2 the non-secular operator 
3/8P in (6) can be neglected in the first approximation. Then, the internal 
evolution of the system will be described by the Hamiltonian — \H! d , while 
the evolution in a rotating frame with no alternating field is described by the 
Hamiltonian H! d . 

The transfer to a tilted rotating frame is a formal operation and cannot 
actually influence the system evolution. At the same time, the influence of 
alternating field 90° -pulse on the system is described by the same opera- 
tor that describes the transfer to a tilted frame. Thus, if we combine the 
sufficiently long application of a strong alternating field to the system with 
the application of short pulse, we can create a situation where the internal 
interactions sign in the system Hamiltonian will be changed with an accuracy 
determined by the correlation between u\ and ujl 

1.2 Magic echo 

The density matrix 

a = 1 - (3oj I x (7) 

describes the state of the of the spin system in a strong external constant 
field after the 90°-pulse was applied to the system. Here (3~ l is the system 
temperature. 

The transverse component of magnetization I x (t) causes the free induc- 
tion signal. 

The exciting experiments of Rhim, Pines and Waugh ^ have shown, 
that the free induction signal restores in time which is much longer than the 
spin-spin relaxation time T2 (T2 is defined as the time required for the free 
induction signal decay under normal conditions). Before the experiments ^ 
this time was considered as a thermodynamic relaxation time in spin systems. 
The phenomenon observed in ^ was called "magic echo". 
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The transverse component of magnetization operator I x (t), to which the 
free induction signal corresponds, is non-diagonal in the energy representa- 
tion. Consequently, the 2nd law of thermodynamics should not apply to the 
I x (t) evolution. Meanwhile, if the irreversibility of macrosystems evolution 
does exist, and if we define the time of equilibrium state establishing in the 
system as the time of irreversible disappearance of the non-diagonal matrix 
elements in the density matrix, we require the condition 

?£ = -i[H,v]=0 (8) 

to be fulfilled in equilibrium, which means that the density matrix should 
become diagonal in energy representation. The disappearance of the non- 
diagonal matrix elements in the density matrix at the irreversible process 
of establishing the equilibrium in the system should, logically, also be irre- 
versible. 

The fact that the I x (t) evolution was reversible in the experiment < - 2 - ) made 
the authors come to a conclusion, that the spin temperature concept should 
be treated cautiously. 

Thus, on the one hand, the spin temperature concept is confirmed by a 
huge number of experimental data, and it is hard to doubt that the spin 
system energy levels are populated according to the Boltzmann distribution. 
On the other hand, the non-diagonal operator I x (t) in the system density ma- 
trix does not irreversibly disappear during a time period which considerably 
surpasses T 2 . If we assume that the spin temperature concept is correct, and 
that the spin system evolution is really irreversible, then the experiments ( 2 ) 
can lead us to a conclusion that T 2 is not a thermodynamic relaxation time. 

If the I x (t) evolution in experiments ^ turned out to be irreversible, 
the authors of ^ would not have had a basis for claiming that the spin 
temperature concept should be treated cautiously. But on the other hand 
the I x (t) evolution in that case would not have been describable on the basis 
of mechanics approach. Such result would have directly demonstrated the 
incompleteness of the quantum theory and the significance of it would have 
been equivalent to this fact. 

The operator I x = Ixi has a simple structure and is a sum of one- 

i 

particle operators. The modern pulse NMR has the methods ^ which allow 
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to bring the spin system to a state described by the density matrix 

a = l-pH' i (9) 

with high value of the inverse temperature f3. 

When # y -pulse is applied to the spin system, the density matrix is trans- 
formed to 

a R = 1- /3[]-{3 cos 2 6 -l)H' d + lp sin 2 6 -^-Q sinO cos6] (10) 
2 8 4 

where 

Q = S ^jVM+j + I-j) + hj{I+i + I-i)]- 

i<j 

and P had been determined by (5). 

Operators P and Q correspond to the interactions of all particles of the 
spin system. It would be natural to expect that the evolution of those op- 
erators would be different from the operator I x evolution. The paper ® 
studied the behaviour of the operator P and Q under the condition of the 
time sign change. In energy representation, Q has only non-diagonal matrix 
elements, but the energy reservoir corresponds to the operator (3/8) P at the 
time of alternating field being applied. The experiments were carried out 
on the nuclei spin system of 19 F in a single crystal CaF 2 . A dipole magic 
echo was discovered - the restoration of the signal confined to the dipole 
interactions in the time that considarably surpasses T%. It turned out that 
the peculiarities in the behaviour of the operator (3/8)P to which the energy 
reservoir corresponds cannot be explained on the basis of reversible equations 
of mechanics. 

The purpose of paper ( 3 ) was to find the dipole magic echo and to com- 
pare the behaviours of the magic echo signal from P and Q operators. The 
measurements were performed in one crystal orientation with respect to the 
external constant magnetic field at a relatively small range of the alternating 
field amplitudes. 

Very important conclusions regarding the irreversibility can be made from 
the results of ®. That is why the need for the continuation of the time 
reversing experiments in spin systems is obvious. Besides, some aspects of 
^ require correction. 
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2 Experiments, discussion 



In this paper we observe and study the dipole magic echo in a wide range 
of external alternating field values at various crystal orientations relative to 
the constant magnetic field. 

The measurements are performed on the CaF2 single crystal, which is 
very convenient for this type of experiments. The avarage value of dipole- 
dipole interactions in different orientations is proportional to local field ojlIi-, 
where ujl is determined by ojl = (M 2 /3) 1 ^ 2 Using the value of the sec- 
ond moment M 2 of the NMR line for CaF 2 , given by Abragam ( 4 ) we find 
in the case of [111], [110] and [100] orientations uj l /i= 0.88, 1.25, 2.01 G 
respectively. 

The pulse sequences which is similar to those presented in ( 3 ) were used 
to achieve the goals. The temperature was 300 K. The time of nuclei spin- 
lattice relaxation was 8 s. Since the duration of the time reversing part of 
the applied pulse sequences was not more than 1 ms, so the analized system 
may be considered isolated from the lattice during this time. 

2.1 Evolution of the (3/8)P dipole subsystem under 
the conditions of time reversion 

We studied the operator P evolution under the conditions of the time rever- 
sion by pulse sequence 1 (Fig.l). 
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The alternating field phase was changed to opposite one in the middle of 
the time interval of the application of this field. The adiobatic demagneti- 
zation in a rotating frame brings the system under consideration to a state 
which is described by the density matrix (9). The system Hamiltonian in a 
tilted rotating frame during the time of alternating field application is (5), 
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and the initial density matrix is 

°(0) = i- p(Ip-\k) (li) 

If the spin system obey the thermodynamic laws, it should be expected 
that after the equilibrium is established in the system the density matrix will 
look: 

a eq = l-(h{o> 1 I x + lp-±H' i ) (12) 

The experiments in laboratory ^ and a rotating ( 5 ) frame have demon- 
strated, that the unified spin system temperature is established in two stages. 
At the first stage, the Zeeman reservoir and the non-secular dipole-dipole in- 
teractions reservoir (to which the operator (3/8)P corresponds in our case) 
undergo the quick warm mixing. At the second stage, the temperature of 
the newly created subsystem and of the reservoir of the secular part of the 
dipole-dipole interactions get levelled at a slower pace. 

2.1.1 Description of the operator P evolution by reversible equa- 
tions 

Let's review the spin system evolution at the application of the pulse sequence 
1 on the basis of expression (2). The transformation which corresponds to the 
90°-pulse effect and to the transfer to the tilted rotating frame compensate 
each other. That's why by the time the 45°-pulse is applied we get: 

°(\ti) = ^(o)^r 1 

A x = exp (-iH'J±)A-A + 

A ± = expi-i^I, + |p - \K)\} (13) 

If the 45°-pulse is excluded from the pulse sequence 1, then after the 
alternating field is removed the spin system signal is absent. The signal 
observed after the 45°-pulse is determined by the transverse component I y , 
and we can write (I y ) = (I y )i + (I y )2- The value {I y }2 is the contribution to 
(I y ) which is bound to the density matrix operator —\Ti' d . The results of ( 5 ) 
demonstrate, that this contribution is constant under the conditions of our 
experiments and can be easily subtracted from the signal under observation. 
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If u\ ^> ul, then in the first approximation we can write 

P(h +t) = A 2 PA 2 l 

A 2 = exp{-iH' d (t-h 1 )} (14) 
It is demonstrated in ®, that in the case of ti = from (14) it follows: 
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(I y )i = ^/3Tr(li)-G(t) (15) 

where G(t) determines the form of the free induction signal. 

A signal whose amplitude does not depend upon the alternating field 
application time corresponds to (15). 

Fig. 2 presents the dependence of the amplitude of the signal correspond- 
ing to the (3/8)P operator in the density matrix upon the alternating field 
application time t±. 

It follows from the Fig. 2 that the signal in the pulse sequence 1 decays as 
t\ grows. The signal decay turned out to grow with the transition from the 
[111] to [100] orientation, but not depend upon the iO\ in every orientation. 

If t\ = N-^, then, using the Magnus expansion ^, we find: 

3 1 

exp{-i(cji/ z + -P - 7:^d)0 = exp(-iu 1 I z t) exp(-iFt), (16) 

o Z 

F = + V 1 Uk 

The Tik operators have the dipole-dipole interactions magnitude in the 
k+1 power. The correction to the average Hamiltonian — \Ti' d , corresponding 
to the first term of the sum by k in (16), is 

W (D = ( 3 )2 [^~ 2) ,^ 2) ] = 3 1 [n' d ,n d - 2) - nf] 



2oo x "82 2uo x 

We shall not need an explicit form of the corrections to — \H! d of a higher 
order. 

If the alternating field phase changes often, like in ( 2 \ the odd power cor- 
rections to —\H'd in (16) vanish ( - 2 ' 6 - ) . In our experiments, the correction TiP^ 
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does not disappear, but the signal decay contribution due to decreases 
because of the phase change of alternating field. 



<Iy>\m 



0.75 - 



0.50 



0.25 - 




250 300 350 400 450 500 
fi, (/is) 

Fig. 2. Dependence on t\ of the amplitude of the magic echo signal due 
to the operator P at ^ = 12.5G(D); ^ = 25.3G(o); ^ = 57.2G(A) in the 
[100], [110] and [111] orientations 



Using (16), by the time the 45°-pulse is applied we get (in the case of 
constant alternating field phase): 

P{\tx) = MP A? (18) 

A 3 = TexpH f 1 exp(~H' d t)n? ■ exp^H' d t)dt] 
JO 2 2 

Writing down (18) we neglected the terms with k > 1 in the expression 
for F (16), because under the conditions of our experiments the ratio of 
— \7~L'd 1° the (fc+i)i (2^! 1 fc ) fc * s °f ^ e or der °f 10 4 already for k = 2. Hence, 
the correction terms with k > 1 cannot contribute considerably to the decay 
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of the signal which is due to the operator (3/8) P. We ommitted also the 
non-secular operator Ti.^ in expression for Ti^. 

Taking into account the change of the alternating field phase we have: 

P^h) = AJ>A? (19) 

A, = exp(-iWi|) exp[*(^ + <H?)\] exp[*(^ - <H?)\]. 

2.1.2 Analysis of the alternating field ungomogeneity influence on 
the experiments results 

It was shown experimentally in ( 3 ), that the signal in the pulse sequence 1 de- 
cays much faster in the case of the constant alternating field phase. It would 
be natural to give the following explanation: when the alternating field phase 
is constant, the contribution to the decay grows by means of the operator 

\ Besides, the changing of the alternating field phase reverses the direc- 
tion of the isochromates precession in the rotating frame and compensates 
the field inhomogeneity influence on the signal under observation. 

The time of the alternating field application in both phases was divisible 
by Fig. 2 demonstrates, that the signal decay at pulse sequence 1 does 
not depend upon the alternating field amplitude. At the same time, the 
signal decay time turned out to vary at different crystal orientation. 

If it was operator Ti^ playing the major role in the signal decay in 
pulse sequence 1, then the signal would have been decaying at a slower pace 
when the uo\ would grow. Also, if it was the not totally compensated field 
inhomogeneity that determined the signal dependence on ti, the signal in 
pulse sequence 1 would have been decreasing when uj\ would grow. Since the 
field inhomogeneity contribution to the signal decay is directly proportional 
to u)i, and the operator contribution is in reverse proportion to Ui, it is 
also not possible to explain the signal independence from lo\ by sum of those 
two contributions. 

We have also achieved purely experimental evidence of the fact that the 
signal decay in pulse sequence 1 cannot be connected to the alternating field 
inhomogeneity. Indeed, the alternating field inhomogeneity does not depend 
upon the crystal orientation, and the fact that the signal heavily depends 
on the orientation points out that the signal decay is not determined by 
the field inhomogeneity even at its maximum. Finally, we have taken the 
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measurements where we used the coils of varying size to create the alternating 
field, which, obviously, produced the fields of different homogeneity However, 
the measurement results did not change when the coils were changed. It 
would also be mentioned that the measurements were performed in this work 
and in ( 3 ) using different equipment and different samples. The results for 
the same orientations and with the same value coinsided, which is the 
evidence of there correctness. 

2.1.3 Irreversibility of the (3/8)P subsystem evolution 

The pulse sequence "c" used in ( 3 ) allow to measure the decay of (3/8)P 
operator signal, which is due to the operator — \Ti! d in the Hamiltonian (5). 
This decay time occurs in orientation [111] was equal to 120 /is. 

The operator 7i^\ which, if the quantum mechanics description is cor- 
rect, determines the dipole magic echo signal decay, is on the two orders of 
magnitude smaller than the operator —\H' d , when the values of cu 1 are those 
used in our experiments. Correspondingly, even not taking the alternating 
field phase change into account, the signal decay time t d in pulse sequence 
1 should be of two orders longer than 120 /xs. The alternating field phase 
change reverses the sign of uj\ and this increases the expected time t d even 
more (see (18) and (19)). Thus, there is ground to think, that the decay 
times t d which are observed at pulse sequence 1 and do not exceed 350 /zs are 
much shorter than the decay time that follows from the system description 
based on expressions (18) and (19) which follows from (2). 

Next, at the transfer from one orientation to another, the value of the 
operator changes in proportion to the second power of the local field, i.e. 
in proportion to the change of the second moment M 2 , which grows 5 times 
when the system transfered from orientation [111] to orientation [100] ^ . 
Hence, the difference between the value of the operator Ti^ in the orientation 
[100] when ^ = 12.5G and in the orientation [111] when ^ = 52.6G is 20 
times, and the signal disappearence time should also be on 20 times different. 
Fig. 2 demonstrates, though, that the times under comparison for the signal 
in pulse sequence 1 are not more than 2 times different. 

Finally, the signal in the pulse sequence 1 does not depend on uj\ contrary 
to the predictions of the theory based on equation (2). 

Thus, the behaviour of the dipole-dipole interaction subsystem, to which 
the operator (3/8) P corresponds, under the time reversion conditions cannot 
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be described on the basis of expressions (1) and (2). This fact leads us to 
the following coupled conclusions: 

a) the evolution of the subsystem corresponding to the operator (3/8)P in 
density matrix is irreversible; 

b) the process of a system coming to an equilibrium state really is irreversible 
and cannot be described by the reversible expressions (1) and (2); 

c) the macrosystem under analysis being isolated does not lead to the re- 
versibility of its evolution. 

At the same time, if the subsystem (3/8)P evolution were reversible, 
it would been the evidence either of the fact, that the spin temperature 
concept is inapplicable, or of the fact that the system's thermodynamic re- 
laxation time is much longer than the time of experiment, or of the fact 
that the mechanics equations are fantastically exact and the thermodynam- 
ics irreversibility is illusionary. But the 2nd law of thermodynamics is the 
generalization of the experience and the discussion of its correctness makes 
no sence, while the exactness of system evolution description on the basis of 
mechanics equations is limited by the extent of the mechanical theory com- 
pleteness. The experimental results which we recieved studying the evolution 
of the subsystem of non-secular dipole-dipole interactions, to which the op- 
erator (3/8)P corresponds, cannot be described by the reversible mechanics 
equations and, hence, demonstrate the incompleteness of those equations. 

2.1.4 Thermodynamical description of the evolution of the (3/8)P 
subsystem 

We select the irreversible component of the evolution of the (3/8)P sub- 
system using pulse sequence 1. It makes sence to introduce its temper- 
ature j3~ l {t) for the thermodynamical description of this subsystem irre- 
versible evolution. The following integro-differential equation for the in- 
verse temperature (3(t) was obtained in ^ on the basis of non-equilibrium 
thermodynamics methods: 




(20) 



where 



Gi(f'-f) 
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It can be seen from (20), that the rate of the reverse temperature change 
does not depend on ui, which corresponds to the fact that the signal observed 
in pulse sequence 1 does not depend on cui. The G\{t) function cannot be 
calculated explicity. Let's write it down as 

G 1 (t' -t) = (nuj L ) 2 G 2 (t' -t) (21) 

where the value of n is of the order of 1. Let us use for G 2 (£) the Gaussian 
approximation, which is natural for regular magnetic of the CaF 2 type ^: 

G 2 (t'-t) =exp(--M(t' -t) 2 ). (22) 

The value M must be comparable to the second moment M 2 of the NMR 
line in CaF 2 . We get M 2 = 2.55 -10 10 , 0.99 -10 10 and 0.5 -10 10 s~ 2 for the [100], 
[110] and [111] orientations ( - 4 * ) . Since operator —\H! d is in the function G\(t) 
exponent, let's write M = M 2 /4. The dotted lines in the Fig.2 correspond 
to the equation (20) solution when n = 0.45 and when M = M 2 /4 for every 
orientation. The agreement with the experiment is the best when the starting 
point is moved 80 /xs in t in the equation (20). It can be explained by the 
fact, that when time t\ is less than 80 fis, the system stochastization after 
the alternating field application does not show yet. Anyhow, there was no 
considerable signal decay at t\ less than 80 /is. 

A total coincidence of theory and experiment when the relaxation pro- 
cesses are described by the non-equilibrium thermodynamics methods seems 
somewhat suspicious. That is why we decided not to try and select such 
version of the correlation function Gi(t) which would make the calculated 
curves to reproduce the experimental data exactly. Besides, the change in 
the alternating field phase can influence the subsystem (3/8) P thermody- 
namics evolution and introduce changes in the equation (20) description of 
this evolution, especially at small values of t±. 

It can be seen from Fig. 2, that in the case of Gaussian approximation of 
G 2 (t) function, when the parameter values correspond to the sample under 
observation, the eguation (20) solution describe the character of the signal 
dependence on orientation and time very well. The degree of the quanti- 
tative correspondence of the theory ^ to the measurement results may be 
considered quite enough. 
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2.2 Evolution of non-diagonal operator Q under the 
condition of time reversion 

Let us now review the spin system evolution at pulse sequence 2 (Fig.l). 
After the 45°-pulse the system state is described by the density matrix 

° = 1 - \P(K + \ P ~ \Q) (23) 
The operator Q in (23) causes the signal ^ : 

(Iy) = 1^(0 J^it) ( 24 ) 
Similar to (19) we obtain: 

Q{\h) = A 4 QA 4 \ (25) 

In accordance with (25), when t = |ti, a signal confined with the operator 
Q appears without any additional influence on the system. The amplitude 
of this signal turned out to depend upon both the orientation of the sample 
and the alternating field amplitude. Hence, the measurement results are 
presented on Fig.3 - Fig.5. Figs. 3 and 4 demonstrate the dependence of the 
dipole magic echo signal upon the alternating field application time in the 
orientation [100] and [110] when the ratio ^ was changed in the interval 12.5 
- 52.7 G. 
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fi, (/is) 

Fig. 3. Dependence on t\ of the amplitude of the magic echo signal due to 
the operator Q at f = 12.5G(D); <f = 25.3G(o); ^ = 52.7G(A) in the[100] 
orientation 
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Fig. 4. Dpendence on t\ of the amplitude of the magic echo signal due to 
the operator Q at f = 12.5G(D); <f = 25.3G(o); ^ = 52.7G(A) in the [110] 
orientation 
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<Iy> 




100 200 300 400 500 600 

ti, (fis) 

Fig. 5. Dependence on t\ of the amplitude of the magic echo signal due 
to the operator Q at ^ = 52. 7G in the [100] - (□); [110] - (o); [111] - (A) 
orientation 



Fig. 5 demonstrates the signal dependence on t\ with — = 52.7 G at 
various crystal orientations. The measurement results for orientation [111] 
at ^= 12.5 G and ^ = 25.3 G coinside with those in ( 3 ) and thus are not 

7 7 

presented here. 

The facts that the signal in pulse sequence 2 grows when uj\ grows and 
that the signal depends on orientation demonstrate that the alternating field 
inhomogeneity influence on the effects under observation is neglectable small. 
If we consider that it is the operator Hi that causes the decay of the signal 
in pulse sequence 2, then the amplitude should grow when P grows, and 
decrease when uj 1 grows. Figs. 3-5 show, that the experimental data qual- 
itatively correspond to this conclusion. The inference that the operator Q 
evolution is reversible when a time-reversing pulse sequence is applied to the 
system was made in ( 3 \ where the measurements were taken only in one ori- 
entation, based on the signal dependence on oj\. However, in present work 
we have analized the whole set of experimental data that we obtained when 
studying the Q magic echo in various orientation and with a wide range of 
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loi values, and this analysis questions the correctness of above conclusion. 

Indeed, the signal decay time in pulse sequence 2, similar to the pulse 
sequence 1, proves to be much less than that corresponding to presence of 
the operator Ti^ in the expression for A 4 . Besides, the decay times in 
orientation [111] at ^ = 52.7 G and in orientation [100] at ^= 12.5 G are 

not more than 3 times different, while the value of the operator Ti^ on 20 
times differs. Thus, it is difficult to reconsile the idea proposed in ® of the 
operator Q evolution being reversible at pulse sequence 2 with the results of 
our experiments. 

The signal decay time t d in pulse sequence 2 is longer than that in pulse 
sequence 1, but values of td remain in the same order of magnitude. Accord- 
ing to the basic statements of statistical physics, the evolution results for the 
diagonal terms of the density matrix and the non-diagonal ones are totally 
different. Hence, the characteristics of the P and Q operators evolution in 
the equilibrium establishing process should also be different. In our exper- 
iments, this difference appears in the fact that the signal corresponding to 
the operator Q depends on oj\. 

2.3 Comparison of the peculiarities of free induction 
signal magic echo and dipole magic echo 

In ref.*- 2 -*, the alternating field phase was changed every -seconds when the 
time reversion situation was created. As a result, the corrections to the mean 
Hamiltonian introduced by the odd powers of uj\ become a zero, and the 
alternating field inhomogeneity gets compensated well. The noticable free 
induction signal magic echo was observed up to ti = 650 fis. The authors in 
^ see the reason for the signal decay in the influence of the alternating field 
phase-changing pulses non-idealities. Actually, the alternating field ampli- 
tude in ( 2 ) was 100 G, and 500 pulses correspond to 650 fis of the alternating 
field application. At that pulse rate there really is ground to assume that 
the magic echo signal decay is explained by the pulse non-idealities. 

We can not say how a very frequent alternating field phase change influ- 
ences the irreversible component of the system evolution. But the possible 
demonstrations of the system evolution irreversibility under the physical con- 
ditions which are created by the phase changing field, can be different from 
the case of the long alternating field application without a phase change. It is 
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therefore possible, that the time of the feasible irreversibility demonstration 
when a phase-changing field influences the system exeeds the signal decay 
time as a result of the pulse non-idealities. 

In our experiments the pulse non-idealities could not have influenced the 
observed signals considerably. But the time of the dipole magic echo signal 
decay turned out to be of the same order of magnitude as the time of the 
free induction signal magic echo decay. It gives us ground to assume, that 
the evolution irreversibility of the transverse component of magnetization I x 
can contribute to the signal decay in the experiment ( 2 \ 

On the other hand, in our experiments the energy redistribution occurs 
in spin macrosystem both in pulse sequence 1 and pulse sequence 2. This 
redistribution, in accordance with 2nd law of thermodynamics, is irreversible 
and the irreversible change of density matrix diagonal terms corresponds to 
this redistribution. 

The operators, giving diagonal and non-diagonal density matrix terms, 
are expressed by means of the same one-particle operators of nuclear moment 
components. Correspondingly, the irreversibility of the evolution of the den- 
sity matrix diagonal terms may lead to the irreversibility of the non-diagonal 
terms evolution also. 

At the same time, in conditions of the experiments ( - 2 - ) , the energy redistri- 
bution does not occur in spin macrosystem. As a result, the character of the 
operator Q and I x evolution in time reversion experiments may be different. 

We would like also to stress here, that operators P and Q are many- 
particle ones and their evolution may be distinguished essentially from the 
evolution of operator I x because of this reason. 

We see that the comparison of the paper ^ and our experiments results 
give rise to many questions. The answers on this questions may be given by 
new experiments only. 

3 Conclusion 

The study of the dipole magic echo that we have performed in a wide range 
of the values of lo\ and at various crystal orientations, allow us to make the 
following main conclusions: 

1) It turned out to be impossible to describe the evolution of the (3/8)P 
subsystem in the time reversion situation on the basis of expressions (1) and 
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(2). The process of the system transfer to the equilibrium state, which is de- 
scribed by the density matrix (12) in the tilted rotating frame, is irreversible. 

2) The dependence of the dipole magic echo signal confined with the 
non-diagonal operator Q upon cu 1 and upon the value of the dipole-dipole 
interactions corresponds qualitative to formula (25) which is derived from 
(2). But the quantitative evaluations demonstrate, that the signal decay in 
pulse sequence 2 happens much faster than it follows from the description 
based on the reversible expresions (1) and (2). 

Probability assumptions are used in a apparent or hidden form when 
a transfer from the reversible mechanics equations to the irreversible non- 
equilibrium thermodynamic ones takes place. As a role, the introduction of 
those assumptions is explained by the facts, that the mechanics equations 
can not be solved exactly and that it is necessary to turn to the shortened 
description of the non-equilibrium processes in the macroscopic systems. The 
experimental and theoretical research performed in ( 3 \ ^ and in this paper 
demonstrate, that it is not the matter of the impossibility of solving the 
reversible mechanics equations exactly, but the matter of their inapplicability 
to the non-equilibrium processes in the macrosystems which obey the 2nd 
low of thermodynamics. This circumstance will allow to take a new view at 
the problem of chaos, both the classical and quantum ones. 

The research that we have performed by far does not exhaust the unique 
possibilities which the dipole magic echo phenomenon presents for the study 
of the correlation between the reversibility and irreversibility in the macrosys- 
tems evolution. We believe that the continuation of the time reversing ex- 
periments in the spin systems will bring a new and, possibly, unexpected 
results. 
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